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I. INTRODUCTION 

The Witten proposal that symmetric deconfined u,d,s- 
quark matter may be the absolute ground state of mat- 
ter [l| has aroused much interest, and the properties of 
strange stars have been widely studied since then. An im- 
portant question is whether the observed compact stars 
are neutron stars or strange stars, which are made up 
of deconfined u,d,s quarks. With the launching of the 
new generation of X-ray detectors Chandra and XMM, 
it is becoming possible now to have accurate measure- 
ment of the radii and surface temperature of compact 
stars 0, 01 El However, theoretical calculations us- 
ing the MIT Bag Model equation of state (EOS) show 
that the mass and size of a strange star are comparable 
to those of a neutron star [(J . Hence it is important 
to identify other observables that can be used to distin- 
guish a strange star from neutron star. One possibility 
is to study the cooling of an isolated compact star. The 
cooling curves of quark matter and neutron matter are 
found to be significantly different due to the difference 
in their thermal properties and energy loss mechanisms 

In this paper, we study the effects of the QCD phase 
transition on the cooling of a compact star and possible 
signatures of the quark phase. Regardless of the valid- 
ity of Witten's proposal, the formation of quark-gluon 
plasma should be favoured in high temperature and den- 
sity we therefore suggest that a strange quark star 
may be formed just after a supernova explosion, in which 
both conditions may be satisfied jl3| . Because the initial 
temperature is so high [l4| Tj ~ 40 MeV, the initial com- 
pact star is likely to be a bare strange star ^(j • When it 
cools down to the phase transition temperature T p , the 
quark matter may become energetically unstable com- 
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pared to nuclear matter, and the strange star will convert 
to a neutron star. During the phase transition, a large 
amount of latent heat, of the order 10 53 erg, can be re- 
leased, which can be a possible energy source of Gamma 
Ray Bursts (GRB's). 

The latest lattice QCD calculations of T p [H [l^l in- 
dicate, though with relatively large uncertainties at high 
chemical potential, that T p drops from its zero density 
value of 140 MeV to about 50 MeV at 1.5 times nuclear 
matter density po = 0.17 fm~ 3 and down to a few MeV 
for density a few times pq. Some previous proto-neutron 
star evolution calculations indeed show that it is feasible 
to reach the phase transition in supernovae ^tJ- While 
there are still large uncertainties in both high density 
QCD and the proto-neutron star evolution, we believe 
it is worthwhile to study the possible consequences of 
the QCD phase transition in supernovae. We assume a 
constant T p in density throughout the star and present 
results for T p = 1, 10 MeV for comparison. We adopt the 
simple picture that matter at temperature above (below) 
T p is in the quark (hadronic) phase. 

It has been argued that the rapid cooling of strange 
stars by pion and e~e + pair production can be a power 
source of GRB's H H 111 E3- However, if the QCD 
phase transition is not considered, the temperature of 
the compact star drops rapidly and the duration of the 
burst - being less than 10~ 2 s - is too short to account for 
long duration GRB's. In our model, the star stays at the 
phase transition temperature for a relatively long time, 
and the photon luminosity is maintained in the range 
10 48 - 10 54 erg s _1 with duration 10~ 3 s up to 10 4 s (or 
even longer) , which is consistent with both long and short 
GRB's. Thus the latent heat may solve the problem of 
the energy supply of GRB's. 

Another signature of quark stars that has been dis- 
cussed is the relatively low luminosity of 7-rays when 
the quark star cools down Our model shows that 

the luminosity in the 7-ray range can be as high as 
10 38-39 erg s with a very long duration of ~ 10 13 ~ 14 
s, which is much longer than the case considered by Page 
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and Usov The satellite Integral launched recently 
is just sensitive to 7-rays in the energy range of interest. 
If a long duration 7-ray source is detected, it could be a 
signature of the phase transition [l8L fli^ . 

A remarkable feature of a strange star to neutron star 
phase transition is the emission of a second neutrino burst 
after the supernova explosion. A similar scenario was 
also proposed in Aguilera et a/.'s work |20j, but with a 
different physical mechanism. In our model, the burst is 
due to the phase transition from a quark star to a neutron 
star, which has a higher neutrino emissivity. In Aguilera 
et aZ.'s work, the burst is due to the trapping of neutrinos 
when the temperature is high, which are then released 
suddenly when the quark star cools. Nevertheless, both 
works suggest that a second neutrino burst is a signature 
of the existence of quark stars. 

This paper is organised as follows. Section II describes 
the EOS used in the quark phase. The stability of strange 
quark star is investigated in Section III. The effect of 
strange quark mass is studied in Section IV. We describe 
how to calculate the cooling history and phase transition 
of compact stars in Section V and VI. Section VII de- 
scribes the models we used. The calculated results for 
the various models are presented in Section VIII. Section 
IX is a short discussion and summary of our work. 



matter properties in the high density regime. Moreover, 
it is worthwhile to study the dependence of compact star 
properties on models of EOS. 

The chemical equilibrium of strange quark matter is 
maintained by the weak-interaction reactions: 



(1) 



and 



s <-> u + e + v e 



s + u «-► u + d- 



(2) 



(3) 



Given the thermodynamic potential of each species 
f2j(i = u, d,s, e), the number densities can be obtained 
from the thermodynamic relation: 



(4) 



where /i^ is the chemical potential. The conditions of 
chemical equilibrium are: 



(5) 



II. COLD EQUATION OF STATE FROM 
PERTURB ATIVE QCD 

Various EOS's have been used to study the proper- 
ties of strange stars. The most widely used one is the 
MIT Bag Model due to its simple analytic form |g. It 
has been pointed out that it is difficult to distinguish the 
strange stars described by the MIT Bag Model from neu- 
tron stars due to the similiarities in their maximum mass 
~ 2Mq (Mq is a solar mass) and radius ~ 10 km [2l| . 
On the other hand, Fraga et al. |21| studied quark star 
structure by using the EOS of perturbative QCD for cold, 
dense quark matter up to order a 2 s , using a modern deter- 
mination of the running of the coupling constant. Their 
results show that strange stars can have a radius of about 
5.8 km and a mass of a typical neutron star, and such 
a small compact star can actually be distinguished from 
neutron stars. Since the strong coupling constant be- 
comes small in the high density limit, perturbative QCD 
may be a fair description of matter in the interior of a 
compact star. Hence we follow Fraga's work and examine 
in details the conditions of absolute and global stability 
of strange quark stars, as well as the effects of strange 
quark mass. It has been argued that perturbative QCD 
may fail to describe the matter at the surface of a strange 
star due to the strong coupling at low density. However, 
the main structure of a compact star is determined by 



u d = Hs- (6) 
Together with the charge neutrality condition: 



2 1, 

o n u - -r (n d + n s ) - n e = 0, 



(7) 



the thermodynamic properties will be determined by one 
independent choice of chemical potential only, which we 
have chosen to be /1 = = fi s in our calculations. The 
total pressure is given by: 



i 

and the total energy density is: 



(8) 



(9) 



At the zero quark mass limit, £l u — fid — fl s implying 
n u = rid — n s . Hence the charge neutrality condition 
is automatically satisfied, without any need of electrons. 
The perturbative QCD thermodynamic potential at zero 
temperature has been calculated up to order a 2 s [23, 123I . 
which in the modified minimal subtraction scheme 1231 
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where G = G -0.5367V/ + TV/ In AT/, G = 10.374±0.13, 
iV/ is the number of quark flavors, A is the renormaliza- 
tion subtraction point, and 
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with u = ln(A 2 /A|^), fa = 11- 2N f /3, fa = 51 - 
19AT//3 and (3 2 = 2857 - 325A^/27. The boundary con- 
dition of a s = 0.3089 at A = 2 GcV for N f = 3 gives 
A^ = 365 MeV. 

All the thermodynamic properties can be obtained 
from the thermodynamic potential if A is fixed. It is 
believed that A//i lies in the range between 2 and 3 [2l| . 
Fig. ^ shows the total pressure of strange quark matter 
relative to the pressure of an ideal gas, Pq, as a function of 
the chemical potential /t. Both the first and second order 
terms decrease the pressure of the strange quark matter 
relative to the ideal gas. It has been pointed out that 
using perturbation theory at zero pressure is invalid |2l| . 
Numerical calculation shows that at zero pressure, a s /ir 
lies between 0.207 and 0.191, which is less than 1 if A//t 
lies between 2 and 3. The EOS's for A//t = 2.473,2.88 
and free gas are shown in Fig. [21 It turns out that this 
EOS is very similar to the MIT Bag Model EOS, with an 
effective Bag constant |2l|, and we would have obtained 
basically the same results using the latter. None of our 
results in the cooling calculation depends on the validity 
of perturbative QCD in the compact star regime. 



baryonic number density at zero pressure is about 1.52no, 
where uq is the normal nuclear matter number density. 
Therefore, strange quark matter can be the true ground 
state, and bare strange stars can exist if A > 2.88/t. 

A weaker condition for strange quark matter to be sta- 
ble is that the average energy per baryon of the bulk mat- 
ter is less than that of 56 Fe, ie., the binding energy per 
baryon of the bulk strange quark matter £binding(SS) is 
larger than that of 56 Fe, £bmdin g (Fe): 



£bii 



5 (SS) 



Mb - Mg 
A 



bindin, 



: (Fe) 



,525 MeV, 
(13) 

where Mb and Mq are the baryonic mass and gravita- 
tional mass of strange quark matter respectively, which 
arc to be calculated by the Tolman-Oppenhcimer-Volkov 
(TOV) equations. We have plotted the binding energy 
per baryon of the maximum mass star against A in the 
lower panel of Fig. |3| We found that the strange star 
can be globally stable compared to iron at infinity if 
A > 2.35//. 

The weak stability condition above overestimates the 
stability of a strange star, since the gravitational binding 
of bulk 56 Fe matter is ignored. A fairer condition is to 
compare the binding energy of a strange star with that of 
a neutron star. There are inevitably some dependences 
on the EOS's used in such a comparison. In Fig. Q] we 
show the parameter regime in which strange stars are 
more stable than neutron stars using one particular EOS 
(HV). We have used several other EOS's to reach a sim- 
ilar conclusion, that strange stars are more stable than 
neutron stars for A > 2.7/1 [Tabic U]. 



III. STABILITY OF STRANGE QUARK 
MATTER 

A strong condition for strange quark matter to be the 
absolute ground state of cold matter is that the energy 
per baryon at zero pressure £q(0) is less than that of 
56 Fe: 

S Q (0) < £ Fe = 930.4 MeV. (12) 

The upper panel of Fig.[3]shows that strange quark mat- 
ter is absolutely stable for A > 2.88/1. At A = 2.88/t, the 



IV. EFFECT OF MASSIVE STRANGE QUARKS 

In the density range of a quark star, i.e. chemical 
potential ~ 300 — 600 MeV, the strange quark mass 
(~ 150 MeV) may alter the EOS. We consider the cor- 
rection of the thermodynamic potential due to strange 
quark mass m s up to first order in a s . The masses of 
u and d quarks are small and can be neglected. The 
individual thermodynamic potentials are Q: 
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where / = u or d, and er is the renormalization point for 
the strange quark mass. It is found that a suitable choice 
for cr is 313 MeV^. 

The EOS's are calculated for A = 2.473, 2.88 and var- 
ious m s numerically [Fig. 03 25]. For m s < 225 MeV, 
the changes in pressure compared to massless EOS are 
less than ~ 5% in the high energy density limit. In the 
low energy density limit, the changes are significant be- 
cause the strange quark mass is not small compared to 
the chemical potentials. Since the global structure of a 
compact star is mainly determined by the high density 
regime - about several times the normal nuclear energy 
density - of the EOS, the effect of the strange quark mass 
on the strange star structure is small. The correction to 
the maximum mass of quark stars due to a quark mass up 
to 150 (225) MeV for each series of A is less than 3(10)% 
[Table [H] . Since the experimental data shows that the 
mass of strange quark is ~ 150 MeV, we can safely ignore 
the quark mass when calculating the global properties of 
a quark star. 

A recent study shows that quark matter in the Color- 
Flavor Locked (CFL) phase can be electrically neutral, 
even though the quark masses are unequal. Electrons 
are not needed to maintain the charge neutrality [26| . 
However, the CFL phase occurs at the very high density 
regime. Electrons or positrons may still be present near 
the surface of a quark star, and they are important for its 
cooling behaviour. More detailed discussion of the effect 
of massive strange quarks and i ts p hysical implications 
can be found in Ref. |U HUlf . 



V. PHASE TRANSITION FROM STRANGE 
STARS TO NEUTRON STARS 

It has long been suggested that strange stars can be 
formed from a phase transition of neutron stars to strange 
stars due to an abrupt increase in density |29l y(j ■ How- 
ever, from the theoretical point of view, formation of 
quark-gluon plasma is favoured when both temperature 
and chemical potential are high enough [L^j. Hence it 
is reasonable to suggest that strange stars are formed 
in supernovae where both high temperature and density 



are achieved [l3j. If the strange quark matter is abso- 
lutely stable for high density and zero temperature, i.e., 
A//i > 2.88, the quark star will remain in the quark 
phase even when it cools down. If the strange star is 
energetically less stable than a neutron star below some 
temperature, it will cool down to the phase transition 
temperature, T p , and change into a neutron star contain- 
ing ordinary baryons. If the baryonic mass Mb is con- 
served during the phase transition, the total conversion 
energy £ conv released is given by: 

E conv = [M G (SS) - M G (NS)]c 2 , (17) 

where Mg(SS) and Afc(NS) are the gravitational masses 
of the strange star and neutron star respectively |3lj| . 
Whether a phase transition can occur and how much en- 
ergy is released depend on both the EOS's of quark mat- 
ter and nuclear matter. We are interested in the possibil- 
ity that strange quark matter is only stable for T > T p , 
and so we choose a A/fi < 2.7, so that when the hot 
strange star cools to low temperature, it will convert to 
a neutron star. For A//i = 2.473, the maximum gravita- 
tional mass is 1.516M© with a baryonic mass of 1.60A/ Q 
and radius 8.54 km. We will use this set of parameters 
in the calculation of the cooling behaviour because the 
maximum mass is close to observational data of compact 
stars. In fact, we have used other values of A//i, and 
the cooling behaviour is qualitatively similar, as long as 
the star undergoes a phase transition. The cold EOS 
of the neutron stars will only determine the conversion 
energy. The cooling behaviour is determined by the cool- 
ing mechanisms and heat capacities in our calculations. 
We simply choose the HFV EOS based on the relativistic 
Hartree-Fock approximation j33| to calculate the conver- 
sion energy. A typical neutron star with Mb = 1.60A/q 
and Mq = 1.40M Q is chosen |32| in the cooling calcula- 
tions. For comparision, the energy released for different 
A//i, together with several commonly used neutron star 
EOS's are summarized in Table ITTT1 For A/jtz = 2.473, 
typically 10 53 erg is released during the process, which 
depends only weakly on the nuclear matter EOS. We 
propose that it can be an energy source of Gamma Ray 
Bursts (GRB's) [25|. 
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VI. COOLING PROPERTIES 

It is believed that in an ultrarelativistic heavy-ion col- 
lision, a hot quark-gluon plasma is formed initially, which 
then cools down to the phase transition temperature T p 
and goes into the mixed phase, if the QCD phase tran- 
sition is first order. When all the quark matter has 
hadronized, the temperature drops again [23] . We borrow 
this idea to describe the cooling of a strange star. When 
the strange star is born , th e temperature can be as high 
as a few times 10 11 K 33]. The surface is so hot that 
all the materials, other than strange-quark matter, are 
evaporated leaving the strange star nearly bare without 
any crust [10| . Since the thermal conductivity of strange 
quark matter is very high and the density profile of the 
strange star is very flat, we will take the uniform tem- 
perature and density approximation. Hence the strange 
star cools down according to the equation: 



a 



dT 
~dt 



(18) 



where C q is the total heat capacity of all the species 
in quark matter, and L q is the total luminosity of the 
star. When the temperature drops to T p , the strange 
star undergoes a phase transition and a latent heat E conv 
is released. We simply take T p to be constant in density 
throughout the star, which has a much smaller value than 
the zero chemical potential value of about 150 MeV be- 
cause of the high chemical potential. During the phase 
transition, we assume that the quark and neutron mat- 
ter are distributed uniformly and calculate the luminos- 
ity of the mixed phase by the weighted average of that of 
the quark matter and the neutron matter [Section IVI Dj . 
When the strange star has lost all its latent heat and con- 
verted completely to a neutron star, it then follows the 
standard cooling of a neutron star with an initial temper- 
ature of T p . The detailed thermal evolution is governed 
by several energy transport equations. We adopt a sim- 
ple model that a neutron star has a uniform temperature 
core with high conductivity and two layers of crust, the 
inner crust and the outer crust, which transport heat not 
as effectively as the core or quark matter. The typical 
thickness of the crust is ~ 10% of the radius, and we 
can use the parallel-plane approximation to describe the 
thermal evolution of the inner crust. The thermal history 
of the inner crust can be described by a heat conduction 
equation: 



Ccrust ' 



dT 
dt 



d_ 

dr 



K 



dT 
dr 



(19) 



where c crus t is the specific heat of the inner crust, K is 
the effective thermal conductivity, and e u is the neutrino 
emissivity. As a rule of thumb, the effective surface tem- 
perature T s and the temperature at the interface of inner 
and outer crust, Tf,, are related by |35| : 



1.288(7^14) 



0.455 



(20) 



where g s \4 is the surface gravity in the unit of 10 cm 
s~ 2 , T{,g is the temperature between the inner and outer 
crusts in the unit of 10 s K, and T s § is the effective surface 
temperature in the unit of 10 6 K. The luminosity at the 
stellar surface, £ S urfacc, is equal to the heat flux at the 
interface of inner and outer crusts: 



K 9T L 
or 



rf'acc 



(21) 



where R is the radius of the star. The boundary condition 
at the interface of the core and inner crust is: 



dT dT 

r i i' a I core 



(22) 



dt dr 

where C corc is the total heat capacity of the core, A COIC is 
the surface area of the core, and L™ rc is the total neutrino 
luminosity of the core. 



A. Heat capacity of quark stars 

The total heat capacity is the sum of the heat capac- 
ities of all species in the star. Without the effect of su- 
perfluidity, the quark matter can be considered as a free 
Fermi gas, and the specific heat of quark matter is given 
by |Mll: 



c q = 2.5 x 10' 



\p/ Po ) 2/3 T 9 erg cm" 3 K" 



(23) 



where p is the baryon density and po = 0.17 fm~ 3 is the 
nuclear matter density. Some authors may use energy 
density instead of number density. The difference is only 
of a numerical factor and does not change the cooling 
curve much. In the superfluid state, the quarks will form 
Cooper pairs. The specific heat will be modified as |38[ 



3.15e Q 



„sf 



2.5 - 1.66T + 3.64T 2 







for 0.2 < T < 1 

for f < 0.2, 

(24) 

where k B T c = A/1.76, f = T/T c and A is the energy 
gap in MeV in BCS theory. It has been argued that for 
quark matter, even with unequal quark masses, in the 
Color-Flavor Locked (CFL) phase in which all the three 
flavors and colors are paired, quark matter is automat- 
ically charge neutral and no electrons are required |26j| . 
However, for sufficiently large strange quark mass and the 
relatively low density regime near the stellar surface, the 
2 color-flavor SuperConductor (2SC) phase is expected 
to be preferred. Therefore in a real strange star, elec- 
trons are believed to be present. The contribution of 
electrons can be parametrized by the electron fraction 
Y e , which depends on the model of strange stars. We 
choose Y e = 0.001 as a typical value. The specific heat 
capacity of electrons in the strange star phase is given by 



YeP 



2/3 



r, = 1.7 11)-" T g erg cm~ 3 K _1 . (25) 
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The heat capacity of electrons is unaffected by the su- 
perfluidity of quark matter. Hence it dominates the total 
heat capacity of the strange star when the temperature 
drops below T c . 



B. Luminosity of quark stars 

The total luminosity is the sum of contributions from 
all energy emission mechanisms, including both photon 
and neutrino emmission. The cooling of quark-gluon 
plasma has been studied for many years [221 an d can 
be divided into fast and slow cooling processes. Two of 
the most popular fast cooling processes are the electron- 
positron pair production and quark URCA process. Two 
well known slow cooling processes are thermal equilib- 
rium and non-equilibrium blackbody radiation. We have 
included all these cooling processes in the our calcu- 
lations. The recently proposed pion production in a 
strange star is found to be a very effective cooling mech- 
anism, and it may explain the energy supply of Gamma 
Ray Bursts (GRB's) Q. Therefore we also studied the 
effect of the pion production cooling. Luminosities of 
various cooling mechanisms are plotted in Fig. and 



1. Pion Production 

It is known that for a hot quark-gluon plasma de- 
scribed by a Cloudy Bag Model, pions can be produced 
through two mechanisms: thermal excitation and col- 
lisions between quarks and the bag surface || 0, l4l| . 
When the pions leave the quark star surface, they will 
decay into photons and e + e~ pairs through the following 
channels: 



27 



v e 



(26) 



(27) 



(28) 



It has been shown that the production of pions is a very 
powerful source of e + e~ pairs and photons, with lumi- 
nosity ~ 10 54 erg s _1 . Such a powerful source of e + e~ 
pairs and photons may serve as the source engine of 7-ray 
bursts. The pion emissivity is estimated to be 



Pt. 



2k B T 



4nR 2 



(29) 



where = 7.1 x 10 31 erg cm -3 is the energy density of 
pion field at the stellar surface, T is the temperature of 
the quark star, which is taken to be uniform, and m n — 
140 MeV is the pion mass. In the superfluid state, the 
collisions between quarks are suppressed due to the pair 
up of quarks, and the pion luminosity is then reduced by 
a factor of exp(— A/T). 



2. Quark URCA Process 

The dominating neutrino emission process is the quark 
URCA process: 



d — > u + e + v ri 



(30) 



u + e — » d + v e . (31) 
The neutrino emissivity was estimated as [3(|: 



e d ~ 2.2 x I0 2b a s 



y e 1/3 T 9 6 erg cm- 3 S - 1 , (32) 



where a s is the strong coupling constant, and we have 
chosen a s = 0.4 as a constant value throughout the quark 
star. Note that the definition of the strong coupling con- 
stant is different from that of Iwamoto's. In the super- 
fluid state, the neutrino emissivity is suppressed by a 
factor of exp(-A/T). 



3. Electron-positron Pair Production 

It has been pointed out that the bare surface of a hot 
strange star is a powerful source of e + e~ pairs due to 
the strong electric field at the surface 0. The e + e~ pair 
production rate is [io|: 



n± ~ 10 3y T 9 3 exp(-11.9/T 9 )J(Os~ 



where 



l£ 3 ln(l + 2£- 1 ) vr 5 
3 (l + 0.074£) 3 + ~6~(13.9 + 4 ' 



a ef 



7r ksT 



0.1 



knT ' 



(33) 



(34) 



(35) 



a = 1/137 is the fine-structure constant and ep = 18 
MeV is the Fermi energy of electrons. The luminosity of 
e + e~ pairs is given by: 



where e± 
particles. 



,F ± ~£ ± n ± , (36) 
^ m e c 2 + ksT is the mean energy of created 

4- Thermal Equilibrium Radiation 



The thermal equilibrium radiation of frequency ui less 
than the plasma frequency lo v ~ 20 — 25 MeV is greatly 
suppressed due to the very high density of the quark- 
gluon plasma[g,u3- The luminosity of thermal equilib- 
rium photons is 10]: 



J eq — 



exp(w/T) - 1 ' 



(37) 
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where 



i r' 2 

= 7T~ / d9 sin 9 cos 9 D(lu, 
2tt Jo 



0), (38) 



D(lu,9) = 1 — (R± + i?||) is the coefficient of radiation 
transmission from the quark-gluon plasma to vacuum, 
with 



R± = 



sin 2 {9-9 ) _ tan 2 (6»-6» ) 
, Jin = 



sm 2 {9 + 9 ) " tan^ 



iQ = arcsm 



sh\9\J 1 — (ujp/uiy 



(39) 



(40) 



5. Non-equilibrium Blackbody Radiation 

The above processes are very powerful sources of en- 
ergy emission and will dominate the cooling process at 
very high temperature. Once the temperature drops, the 
cooling process will be dominated by the relatively low 
power non-equilibrium blackbody radiation 42] : 



L 



ncq 



lQ- 6 L bb , 



(41) 



where L bb — <iirR 2 <jT 4 is the blackbody radiation lumi- 
nosity. 



C. Microphysics of the neutron star cooling 

There are many different models of neutron star cool- 
ing. Since we mainly focus on the cooling of the quark 
phase and examine the phase transition process qualita- 
tively, the cooling of neutron stars can be taken from any 
model available in the literature. We simply adopt the 
model of neutron stars described by Ng 0, |3!j . The heat 
capacities of neutron matter in both normal and super- 
fluid states, and of the electrons are: 

c„ = 2.3 x 10 39 M»p7; 4 2/3 T 9 erg hT 1 for f > 1, (42) 



p + e — > n + v e , 
with the neutrino emissivity 0, 0] : 



(46) 



euRCA = 4.00xl0 27 (y e/9 / jO0 ) 1/3 2t erg cm-V 1 for T > T c , 

(47) 



£ urca = euRCAexp(-A/T) for T <T C , 



(48) 



with Y e = 0.1, electron-proton Coulomb scattering in the 
crust: 

e+p e+p + v + V (49) 

with luminosity pjLl44|: 

L C J = 1.7 x 1 39 M* (M CT /M) Tg erg a -1 , (50) 

where (M cr /M) is the fractional mass of the crust « 5%, 
and the neutrino bremsstrahlung process: 



n + n—^n + n + is + i' 



(51) 



with luminosity jg]: 

L n v n = 4.3 x 10 38 pi4 3 r 9 8 erg s" 1 for T > T c 







for T < T c . (52) 



The surface luminosity will be of the blackbody radia- 
tion with the effective surface temperature T s : 



(53) 



The blackbody radiation will be the dominating cool- 
ing mechanism after the neutrino emission processes are 
switched off. 

For the thermal conductivity of the inner crust, we 
simply choose a temperature dependent model |45j : 

K = 8 " erg cm" 1 s" 1 K . (54) 
Tio 

The choice of thermal conductivity will not be important 
after the epoch of thermal relaxation, which is of the 
order 10 — 100 years. The temperature of the inner crust 
and the core will be uniform after that. 



and 



3.15c„ 
— = — i 
T 



2.5 - 1.66T + 3.64T' 



1.9 x 10 37 M*p^°T 9 erg K 



l/3r 



for 0.2 < T < 
(43) 

(44) 



where is the mass of neutron star in units of solar 
mass, pi4 = p/10 14 g cm" 3 , and we take T c = 3.2 x 10 9 
K for our calculations. 

The neutrino emission mechanisms are the direct 
URCA process: 



p- 



1 D. Handling of the Phase Transition 

The actual situation during the phase transition is very 
complicated, involving detailed hydrodynamical simula- 
tions, how quark matter is transformed into ordinary 
hadronic matter and so on. Many authors have dis- 
cussed the theoretical modeling of the phase transition 
in heavy-ion collisions, neutron star-strange star phase 
transition and the cosmological quark-hadron transition, 
but the uncertainties are of course very large at this stage. 
Without involving the details, we simply assume that the 
phase transition is first order and evolves quasi-statically. 
(45) The quark and neutron matter are distributed uniformly 
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in the mixed phase. We develop a simple model to de- 
scribe the luminosities of the mixed phase in order to 
capture the main features just before and after the phase 
transition. The detailed results during the phase transi- 
tion is of course inaccurate, but the main features should 
be approximately correct, as long as the phase transition 
does not disrupt the neutron star completely. 
The energy loss during the phase transition is: 



The luminosity of thermal equilibrium photons during 
the phase transition is given by Ea. l|37[) where oj p = uj p T . 
We can see that the lower limit becomes smaller dur- 
ing the phase transition. The luminosity reduces to that 
of thermal blackbody radiation at the end of the phase 
transition. Hence we expect an increase of equilibrium 
photon luminosity as a signature of the phase transition. 



-Eloss — 



PT 



PT 



f/f, 



(55) 



where Lp q T and are the total luminosity of the quark 
matter contribution and nuclear matter contribution dur- 
ing the phase transition respectively, tQ T is the starting 
time of the phase transition. We define the fractions 
of total baryonic number of quarks and neutrons in the 
mixed phase to be j4ss and Ans respectively. Ass and 
Ans w iH be related to the conversion energy E conv and 
energy loss during the phase transition £i oss by: 



Ass = 



and 



.4 



NS 



(56) 



(57) 



1. Pion Production, e + e Pair Production 

Since the pion production and e + e~ pair production 
are the features of a hot quark star, we expect the produc- 
tion rates to decrease with the fraction of quark number 
during the phase transition. We simply treat the lumi- 
nosities of both contributions to be: 



2. Thermal Equilibrium Radiation 



(58) 



The plasma frequency is related to the baryon number 
density as [lOj : 



/8tt ■ 



(59) 



~ P ~V3 fx 
where p, ~ Hc^rib) 1 ^ 3 is the chemical potential. Hence 



,V3 



(60) 



For a typical strange star, uj p ~ 20 MeV. Therefore dur- 
ing the phase transition, the plasma frequency will be 
related to the fraction of quark number as: 



l T ~ 20A SS 1/3 MeV. 



(61) 



3. Quark URCA Process 

From Eg. 132(1 . the neutrino emissivity of a quark star 
td oc p. Hence, 

(62) 



PT A 

(d = ( -d,Ass- 



4- Neutrino Emission Process of the Neutron Matter 

For the neutrino emissivity of neutron matter, we use 
the fractional density of neutron matter in the density 
dependent terms of Ea. l|4T|) and JSU), i.e.: 

p -> A NS/ o. (63) 

For the electron-proton Coulomb scattering of Ea. (|50|l . 
we simply calculate its contribution with a weighted 



PT 



NS- 



(64) 

Lattimer et al. have pointed out that direct URCA pro- 
cess can only be switched on in a neutron star at suf- 
ficiently high central density 01 . For a wide range of 
parameters, the required central density is of the order 
of nuclear matter density. We switch on the direct URCA 
process only when the average density of neutron matter 
is above the nuclear matter density. 



5. N on- equilibrium Blackbody Radiation 

Since the non-equilibrium blackbody radiation is not 
the dominating cooling mechanism in the temperature 
range of 10 9 — 10 11 K, it can be neglected in the mixed 
phase. 



VII. THE MODELS 

When a strange star is born just after the stellar 
collapse, its temperature is very high ~ 10 11 K ~ 40 
MeV 14j. We choose an initial temperature Ti = 
40 MeV. The EOS of perturbative QCD with massless 
quarks and temperature correction up to first order in 
a s is |H|: 



P = 



45 



rpi 



2 



1 _ i^f^f. 

4 7T 



-7T 2 T 4 

60 



_ 50 Og 

~ 21 7T 



2 /; 2 



1 

47T 



a, 

1 - 2— 

7T 



(65) 
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Since just after the collapse, T ~ 40 MeV < \x ~ 300 
MeV at the surface, the zero temperature EOS is appro- 
priate for calculating the structure of the star. In our 
cooling model, the only free parameter of the EOS in 
the quark phase is A. It will determine whether a phase 
transition can occur and how much latent heat can be 
released. For a 1.4M compact star, the cooling curves 
are mainly affected by the cooling mechanisms and heat 
capacities of the material, rather than A, which mainly 
affects the structure of the compact star. Therefore we 
only choose one value of this parameter, A = 2.473/1, 
for the cooling calculations. The phase transition tem- 
perature T p involves detailed study of the baryon phase 
diagram, which is still highly uncertain. The latest re- 
sults Il5l indicate that T p drops from its zero density 
value of 140 MeV to a few MeV for density of a few times 
nuclear matter density. Hence we treat T p as a param- 
eter and study the cooling behaviours for different T p . 
We assume a constant T p in density throughout the star 
and present results for T p = 1, 10 MeV for comparision. 
The remaining ingredients are the heat capacities and 
cooling mechanisms of the star. For the quark phase, the 
URCA process, thermal equilibrium and non-equilibrium 
radiation are believed to occur for a bare strange star. 
We include all these processes in our calculations, and 
we choose other combinations of cooling processes to in- 
vestigate the effect of each cooling mechanism. For the 
superfluid phase of quark matter, we choose A = 1 MeV 
for a small gap model and A = 100 MeV for a large gap 
model. The various models presented in this paper are 
summarized in Table llVl 

VIII. RESULTS 

The strange star luminosity L discussed above is typ- 
ically many orders of magnitude higher than 10 37 erg/s 
and may be as high as 10 53 erg/s. In this case, the out- 
flowing wind is optically thick, and at L > 10 42 — 10 43 
erg/s the spectrum of emergent photons is nearly a black- 
body spectrum ^?], |4^. Therefore, for a new born 
strange star, the pion production, equilibrium radiation, 
and electron-positron pair production cannot be distin- 
guished observationally in spite of the fact that the char- 
acteristic gamma-ray energies at these processes near the 
strange star surface differ significantly. They can be dis- 
tinguished observationally only if the phase transition 
temperature is very low, such that the strange star can 
exist in a relatively low temperature. The observables 
are the luminosity and the surface temperature at infin- 
ity, L°° and T s °°, which are related to the stellar surface 
values, L and T s |49| : 

L°°=e 2 ' t '*L, (67) 

where e^ s = yl — 2M/R is the gravitational redshift at 
the stellar surface. The various cooling curves for dif- 



ferent models are shown in Fig. 8.1 - 8.4. Models a to 
c are small superfluid gap models, while Models d to f 
are large superfluid gap models. For small (large) su- 
perfluidity gap models, the phase transition temperature 
considered here is well above (below) the superfluid gap. 
We can get some insight of the effect of superfluidity on 
the signatures of phase transition. 



A. Model a, b 

The cooling curves for Model a and Model b are nearly 
the same. It is because both of them have the same 
thermal properties, and all the cooling mechanisms are 
switched on, except for Model a that the e~e + pair pro- 
duction is not considered. Since for the small superfluid 
gap A = 1 MeV, pion emission is not suppressed for 
T > T c . The pion emissivity dominates the e~e + pair 
production by several orders of magnitude at high tem- 
perature. At low temperature, the cooling process will 
be dominated by the blackbody radiation. Hence both 
models will be dominantly cooled by the same mecha- 
nisms and have the same cooling history. Only if the 
phase transition temperature is very low, we can distin- 
guish them by studying the spectrum, in which Model b 
shows the character of e~e + production. In these mod- 
els, for T p — 1 MeV, the photon luminosity can maintain 
up to 10 50 — 10 54 erg s _1 for about 1 s when the star is 
in the quark and mixed phases. Such a 'long' duration of 
extremely high luminosity, compared to Ng et a/.'s model 
which gives a luminosity that drops below 10 50 erg s _1 
within 0.01 s in the extreme case, is dominated by pion 
emission and maintained by the latent heat of the star. 
This violent fireball easily supplies the energy required 
for GRB's. The neutrino luminosity drops from 10 57 to 
10 48 erg s _1 within 0.001 s and then rises to 10 erg s 
during the mixed phase which lasts for ~ 1 s. 

This scenario of a second burst of neutrinos can be 
compared to two previous similar proposals [171 |2Qj . In 
our model, the burst is due to the phase transition from 
a quark star to a neutron star, which has a higher neu- 
trino emissivity, whereas in previous proposals, the sec- 
ond burst accompanies the phase transition from a neu- 
tron star to a quark star. In Benvenuto's theory, the 
phase transition is delayed by a few seconds after the 
core bounce due to the presence of the neutrinos |17j . In 
Aguilera et al.'s theory [2£|, the burst is due to the initial 
trapping of neutrinos when the temperature is high and 
their sudden release when the quark star cools. If quark 
matter is not as stable as nuclear matter at low tem- 
perature, then there should be yet another phase transi- 
tion back to nuclear matter, which is what we focus on, 
and the "second" neutrino burst we propose is then the 
"third" neutrino burst. 

If multiple neutrino bursts are observed, as may indeed 
be the case for the Kamiokande data for SN1987A [50| . 
whether the compact star changes from the quark phase 
to neutron phase (our model) or the other way around 
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can be distinguished observationally in at least two ways. 
First, our model predicts that the cooling is much faster 
before the phase transition, but it will become slower after 
it. Second, the size of the post-phase-transition compact 
star, being a normal neutron star, would be larger in our 
model. 

For T p — 10 MeV, the high photon luminosity can only 
be maintained for ~ 0.02 s, since the neutrino luminos- 
ity is very high during the phase transition, in which it 
releases the latent heat within a short time. The rises in 
neutrino flux are about two orders of magnitudes, but the 
time scale is so short that they are likely to be masked by 
the diffusion time of neutrinos (~ 1 — 10 s)0| out of the 
dense medium and therefore probably indistinquishable 
from the first burst. 

We have also studied the case of T p = 0.1 MeV (not 
shown in the figures). The photon luminosity can be 
maintained at 10 45 erg s _1 up to 5 x 10 6 s, and it then 
decreases gradually to 10 34 erg s _1 in 10 8-9 s, while the 
neutrino luminosity first drops to 10 34 erg s _1 and then 
rises to 10 41 erg s~ 4 during the mixed phase. If such a 
long duration 7-ray source exists, it will be difficult to 
be explained by other models. For our model with T p < 
0.01 MeV (not shown), the cooling history is similar, but 
with an even longer duration and lower luminosities in 
the mixed phase. 

The cooling curves (of all our models) are basically 
cooling of a bare quark star plus that of a neutron star. 
When the quark star cools to the phase transition tem- 
perature, the cooling curve switches from that of the 
quark star to the neutron star, plus a phase transition 
epoch. We can see that for a typical bare quark star 
without phase transition, the photon luminosity, as well 
as the surface temperature, are higher than that of the 
neutron star in the early epoch. The quark star then 
cools much faster than the neutron star, and it will be 
too cold to be observed. For the cooling of a bare quark 
star with phase transition, there is a significant feature 
different from the case without phase transition. When 
the quark star cools to the phase transition temperature, 
it remains at the phase transition temperature, and it 
releases the latent heat for a long time depending on the 
latent heat and the total luminosity. The total luminosity 
actually depends on the temperature. This means that 
the higher the phase transition temperature, the larger 
the luminosity and hence the shorter the duration of the 
phase transition. The latent heat depends on the choice 
of the EOS's of both the quark and neutron matter, as 
well as the mass of the compact star. Here we have chosen 
A = 2.473^ for quark matter, the HFV EOS for neutron 
matter and a compact star mass of ~ 1.4M©. The latent 
heat is hence 2.08 x 10 53 erg. 



B. Model c 

For Model c, the thermal properties are the same as 
Models a and b. In this model, the pion emissivity is not 



considered, and the dominating cooling process in the 
early epoch is the e~e + production. Hence the photon 
luminosities for Model c in the early epoch is significantly 
lower, and the durations of phase transition for high T p 
models are longer. For T p — 1 MeV, the photon luminos- 
ity drops from 10 53 to 10 48 erg s _1 and can maintain such 
a high luminosity up to 10 s. Such a model seems to be 
in good agreement with long duration GRB's both in the 
energy and time scales. For T p = 10 MeV, the cooling 
curves are similar to those of Model a and b, except for 
the lower photon luminosity with a slight rise near the 
end of the phase transition in Model c. The similarity 
in cooling curves is due to the same dominating cooling 
mechanism of neutrino emission at the high temperature 
epoch in the quark and mixed phase, as well as the same 
cooling behaviour after the phase transition. We have 
again studied the case for T p = 0.1 MeV (not shown in 
the figure) , the photon luminosity can be maintained up 
to 10 38-39 erg s _1 for 10 13-14 s. This model is similar to 
Page and Usov's work [lj , but our work can maintain a 
high luminosity for much longer time. The late and phase 
transition epoch of the cooling curve for T p < 0.01 MeV 
(not shown in the figure) are the same as those of Models 
a and b at the same T p . It is because for temperature 
< 0.01 MeV, the dominating cooling mechanisms in all 
models (including Models d - e as well) in the quark and 
mixed phases are non-equilibrium blackbody radiation. 



C. Model d, e 

As compared to Models a and b, the cooling curves for 
Models d and e are also nearly identical. It is because 
the main effect of a large superfluid gap is to suppress 
pion and neutrino emissivities, which are important for 
the high temperature epoch. Hence the e~e + produc- 
tion dominates the photon luminosity in the early epoch. 
Thus both models are dominantly cooled by the same 
mechanisms. We can see that the effect of superfluidity 
is to make the neutrino bursts more distinct. For T p = 10 
MeV, the neutrino flux rises by over ten orders of magni- 
tudes within a small fraction of second. The effect may 
again be masked by the diffusion of neutrinos. However, 
if T p is as low as 1 MeV, the two bursts of similar flux 
can be separated by as long as 10 5 s, which should be 
observable by modern neutrino observatories. 



D. Model f 

In Model f, e~e + production is not considered. Al- 
though pion emissivity is included, it is again greatly 
suppressed. Hence the cooling history is as if there is 
no fast cooling mechanism. The quark star cools slowly, 
maintaining a relatively low photon luminosity for ~ 10 8 
s for Tp = 1 MeV. One special feature of this model 
is that the photon luminosity rises up for about six or- 
ders of magnitudes at the end of the phase transition 
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for T p = 1 MeV. This is due to the increase in lumi- 
nosity of thermal equilibrium radiation during the phase 
transition [Section TVI Dj . Indeed for Model a - e, if we 
separate the spectrum of thermal equilibrium radiation 
from other photon luminosities, it is also increasing. In 
our models, we expect more photons to be released dur- 
ing the phase transition. The duration is ~ 10 s in this 
model and are shorter for other models. If such a burst 
is observed, it will be a unique feature which is difficult 
to be explained by other models. We remark here that 
detailed radiative transfer calculation may kill or lower 
the peak of this burst due to the surface effect of the 
star. Another special feature of this model is that for all 
Models a - e (except Model c) , the second neutrino burst 
starts within 1 s from the first burst, which may be too 
short to be detected. This agrees with the observational 
data of SN1987a 01 ( we emphasis here that base on the 
one event and the very few neutrinos detected, we re- 
ally cannot tell how many neutrino bursts are there). In 
Model f, the second neutrino burst starts after 10 4-5 sec- 
onds. Such a delayed neutrino burst may be detectable, 
and it may have significant effects on the propagation of 
shock waves in supernova remnants 



1 . Remarks 

In these models, we can see that the size of superfluid 
gap determines how a compact star cools in the early 
epoch. Pion emission is a very efficient cooling mecha- 
nism if the superfluid gap is small. If the superfluid gap is 
large, e~e + production will dominate in the early epoch. 
Also we find that a second neutrino burst is a signature 
of the phase transition, if the neutrino emissivity in the 
neutron phase is higher than that of the quark phase. In 
our models, the URCA process in a neutron star makes 
the second neutrino burst possible. 



IX. DISCUSSION AND CONCLUSION 

We studied the global structure and stability of strange 
stars with the perturbative QCD EOS up to order a 2 s . 
We find that for A > 2.88^, strange quark matter is 
absolutely stable, while for A > 2.35(> 2.7)/x, a strange 
star is globally stable (compared to a neutron star). The 
effect of strange quark mass to the strange star is also 
studied. For a strange quark mass < 150(225) MeV, 
the correction to the maximum mass of strange stars is 
less than 3(10)%. We suggest that a strange star may 
undergo a phase transition to a neutron star when it cools 
down to some temperature T p . 

It has been argued that the rapid cooling of strange 
stars by pion emission can be a power source of GRB [8| . 
However, if the phase transition is not considered, the 
duration of the burst is too short, < 10~ 2 s, which cannot 
explain long duration GRB's, due to the rapid cooling 
without the maintenance of high temperature. However, 



in our model, the latent heat of the phase transition can 
supply the energy of the order 10 53 erg. In our models, 
the photon luminosities can be maintained in 10 48 — 10 54 
erg s _1 with duration from ~ 10~ 3 up to 10 4 s (or even 
longer), which are in the range for both long and short 
7-ray bursts. Thus the latent heat solves the problem of 
the energy supply for the GRB's. 

Another signature for quark stars that has been dis- 
cussed is the relatively low luminosity of 7-rays when the 
quark star cools down For models with low T p , e.g. 
Model c with T p ~ 0.1 MeV, the luminosity due to the 
e~e + pairs in the 7-ray range can be as high as 10 38-39 
erg s _1 with a very long duration of ~ 10 13-14 s, which is 
much longer than the case considered by Page and Usov 
[Tl|. The satellite Integral launched recently is just sen- 
sitive to 7-rays in the energy range of 15 keV - 10 MeV. If 
long duration 7-ray source is really detected, it would be 
a signature of the phase transition. On the other hand, if 
no such source is detected, this does not mean failure of 
our models. Perhaps the phase transition temperature is 
very high, or the latent heat is small. Indeed Integral has 
already discovered a number of soft 7-ray point sources 
in the center of our galaxy |18| . Oaknin and Zhitnitsky 
have also discussed the possibility of the 7-ray sources 
being supermassive very dense droplets (strangelets) of 
dark matter in a recent paper |19| . 

If a second neutrino burst after a supernova explosion 
is detected, it will be a strong evidence for the existence 
of quark star, and in addition, it will support the switch 
on of URCA process, or other fast neutrino emitting pro- 
cesses in neutron stars. However, the absence of the sec- 
ond burst does not kill quark stars, because there could 
simply be no fast neutrino emitting process in the neu- 
tron star phase, or the second burst arises so quickly that 
we cannot distinguish it from the first one. 

It is generally believed that GRB's are related to su- 
pernovae pHl IS . Our model is compatible to the 
supernova connection - a supernova explosion leaves be- 
hind a compact star, which triggers the GRB weeks to 
years later. We suggest that the phase transition from 
a hot strange star to a neutron star may be the central 
engine of GRB's. During the phase transition, the size 
of the compact star increases, and an internal shock may 
develop, which initiates the GRB. 

In our models, the phase transition is treated in an 
oversimplified manner. Detailed study should be made 
which involves hydrodynamics, the EOS, the emission 
properties of matter in the mixed phase and so on. Also 
the duration of the quark phase may be so short that for 
the cooling process, the hydrostatic treatment of quark 
stars may not be appropriate. Hydrodynamic simulation 
may be needed for the whole process starting from the 
supernova explosion. This involves detailed numerical 
treatment which is beyond the scope of this paper. Here 
we discuss semi-quantitatively the signatures left during 
the phase transitions by the presence of the quark phase, 
if it is reached in supernovae. 
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FIG. 1: The total pressure of strange quark matter, P, 
relative to the pressure of an ideal gas, Po, as a function of 
chemical potential, /i. P is calculated up to first and second 
order in strong coupling contant, a s . Here we take A = 2.88/x 
and zero quark mass approximation [21]. 



FIG. 2: Equation of state (pressure, P, vs energy density, e) 
for cold quark matter with massless quark approximation at 
various A, and free gas [21]. 
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FIG. 3: Upper panel: Energy per baryon at zero pressure, 
£q(0), as a function of A//U. The dotted line corresponds to 
the energy per baryon of 56 Fe, which is equal to 930.4 MeV. 
Absolute stability of strange quark matter corresponds to A > 
2.880/4. Lower panel: Binding energy per baryon fbinding(SS) 
as a function of A//x. The dotted line corresponds to the 
binding energy of Fe, which is equal to 8.525 MeV. Global 
stability corresponds to A//i > 2.35/1. 
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FIG. 5: Effect of strange quark mass, m s , on the pressure P. 
Po is the zero quark mass pressure of strange quark matter 
about the same energy density e. 
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FIG. 4: Binding energy per baryon fbinding(SS) as a func- 
tion of A/jj, for baryonic mass of strange stars being equal to 
that of a 1.4M Q neutron star with HV EOS. The dotted line 
corresponds to the binding energy per baryon of the neutron 
star, which is equal to 68 MeV. Global stability compared to 
the neutron star corresponds to A > 2.65/j. 
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FIG. 6: Observed luminosity L°° as a function of temper- 
ature T for small superfluid gap model in the quark phase 
A = 1 MeV. Ltt, L±, L eq , L v and L„ eq denote the pion, e e~ , 
thermal equilibrium radiation, neutrino emission and non- 
equilibrium blackbody radiation luminosity respectively. For 
comparison, the neutrino luminosity of neutron star phase 
is also plotted. SS and NS correspond to strange stars and 
neutron stars respectively. [8] 
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Fig. 8.1 corresponds to Models a and b. Fig. 8.2 corresponds -6 -3 3 6 -6 -3 3 6 9 

to Model c. Fig. 8.3 corrsponds to Model d and e. Fig. 8.4 log (t/ s) log (t/ s) 



corresponds to Model f. In each figure, the left (right) panels 
are for models with small (large) phase transition temper- 
ature T p = 1 (10) MeV. The dotted (dashed) lines are the 
cooling curves of a pure strange (neutron) star without phase 
transition. The solid line represents the scenario with phase 
transition. T a °°, and LJ° denote surface temperature, 
photon luminosity and neutrino luminosity respectively. An 
initial temperature of 40 MeV is assumed. 
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EOS(NS) 


Mb/Mq 


£binding(NS) 

[MeV] 


A/H 


M G (SS)/M e 


£binding(SS) 

[MeV] 


HV 


1.51 


68.4 


2.473 


1.44 


43.5 


HV 


1.51 


68.4 


2.600 


1.41 


62.1 


HV 


1.51 


68.4 


2.880 


1.33 


111.9 


HV 


1.51 


68.4 


3.000 


1.29 


136.8 


HFV 


1.60 


117.4 


2.473 


1.516 


49.3 


HFV 


1.60 


117.4 


2.600 


1.478 


71.6 


rlr V 


i an 

l.DU 


117/1 
1 1 f .4 


z.ooU 


i Ann 


1 1 f .4 


HFV 


1.60 


117.4 


3.000 


1.363 


139.1 


aW+hfv 


1.62 


127.5 


2.473 


/ 
/ 


/ 
/ 


AS BI S F +HFV 

yi BroB +nr v 


1.62 


127.5 


2.600 


1.50 


69.6 


yi BroB + nr V 


1.62 


127.5 


2.880 


1.41 


121.7 


Ji BroB ' AAi v 


1.62 


127.5 


3.000 


1.38 


139.1 


p K24U 
^M78 


1.56 


96.3 


2.473 


1.48 


48.2 


P K240 
^M78 


1.56 


96.3 


2.600 


1.45 


66.2 


r K2A0 
^M78 


1.56 


96.3 


2.880 


1.37 


114.4 


r K240 
^M78 


1.56 


96.3 


3.000 


1.33 


138.4 


G^78 U (NP) 


1.56 


96.3 


2.473 


1.48 


48.2 


gS 2 7 4 8 (np) 


1.56 


96.3 


2.600 


1.45 


66.2 


G^ 2 7 4 8 (NP) 


1.56 


96.3 


2.880 


1.37 


114.4 


G^ 4 8 °(NP) 


1.56 


96.3 


3.000 


1.33 


138.4 



TABLE I: Binding energy, ^binding; f° r various compact stars. A neutron star gravitational mass Mg(NS) = 1.4M Q is assumed, 
and the baryonic masses of strange stars are chosen to be equal those of the neutron stars. The many-body approximation for 
HV, HFV, Ag r B G | F + HFV and 

Gm78° EOS's are relativistic Hartree, relativistic Hartree-Fock, relativistic Brueckner-Hartree- 
Fock + relativistic Hartree-Fock and relativistic Hartree respectively [32] . 



TABLE II: Effect of strange quark mass on the global structure of a strange star. The superscripts of 'max' correspond to 
quantities of maximum mass stars. 







A = 


2.473^ 






m s (MeV) 


max/ N 

e c (eoj 


M max (M ) 


% increase 


i? max (km) 


% increase 





12.0 


1.516 


/ 


8.54 


/ 


75 


12.3 


1.550 


+2.22% 


8.64 


+1.17% 


150 


12.1 


1.533 


+1.11% 


8.59 


+0.59% 


225 


14.0 


1.441 


-4.95% 


8.15 


-4.57% 


300 


16.1 


1.323 


-12.75% 


7.69 


-9.95% 







A = 


2.88^ 






m s (MeV) 


max/ \ 

e c Vo) 


M max (M ) 


% increase 


i? max (km) 


% increase 





8.0 


1.983 


/ 


11.07 


/ 


75 


7.2 


2.022 


+1.93% 


11.33 


+2.35% 


150 


7.3 


1.956 


-1.36% 


11.08 


+0.09% 


225 


8.8 


1.796 


-9.46% 


10.35 


-6.50% 


300 


9.5 


1.641 


-17.26% 


9.92 


-10.39% 
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TABLE III: Total conversion energy Econv released of a compact star for different A. The neutron star gravitational mass 
M G = 1.4M with various EOS's. 



EOS(NS) 


M B /M @ 


a/a* 


M G (SS)/M e 


Sconv/lO^ erg 


HV 


1.51 


2.473 


1.44 


+0.72 


1 1 V 


l.ol 


o Ann 


1.41 


i n i o 
+U.lo 


1 1 V 




9 880 


i-.OO 


— 1.25 


HV 


1.51 


3.000 


1.29 


-1.97 


HFV 


1.60 


2.473 


1.516 


+2.08 


HFV 


1.60 


2.600 


1.478 


+1.40 


HFV 


1.60 


2.880 


1.400 





HFV 


1.60 


3.000 


1.363 


-0.66 


Ag^+HFV 


1.62 


2.473 


/ 


/ 


A« r B l F +HFV 


1.62 


2.600 


1.50 


+ 1.79 


A« r B | F +HFV 


1.62 


2.880 


1.41 


+0.18 


A« BBF +HFV 


1.62 


3.000 


1.38 


-0.36 


U M78 


1.56 


2.473 


1.48 


+1.43 


n K240 
U M78 


1.56 


2.600 


1.45 


+0.90 


r.K240 
^M78 


1.56 


2.880 


1.37 


-0.54 


n K240 
U M78 


1.56 


3.000 


1.33 


-1.25 


G^ 7 4 8 U (NP) 


1.56 


2.473 


1.48 


+1.43 


G^f(NP) 


1.56 


2.600 


1.45 


+0.90 


G»NP) 


1.56 


2.880 


1.37 


-0.54 


G»NP) 


1.56 


3.000 


1.33 


-1.25 



TABLE IV: Models of quark phase presented in this paper. 









Cooling 


Mechanisms 




Superfluid Gap 


Model 


quark 


neq 


eq 


e + e~ 


pion 


A(MeV) 




URCA 












a 


V 


V 


V 


X 


V 


1 


b 


V 


V 


V 


V 


V 


1 


c 


V 


V 


V 


V 


X 


1 


d 


V 


V 


V 


V 


X 


100 


e 


V 


V 


V 


V 


V 


100 


f 


V 


V 


V 


X 


V 


100 



